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Abstract. - We present an exact derivation of the survival probability of a randomly ac- 
celerated particle subject to partial absorption at the origin. We determine the persistence 
exponent and the amplitude associated to the decay of the survival probability at large times. 
For the problem of inelastic reflection at the origin, with coefficient of restitution r, we give 
a new derivation of the condition for inelastic collapse, r < r c = e~ n ' , and determine the 
persistence exponent exactly. 



Consider a randomly accelerated particle, obeying the stochastic equation of motion 

d 2 x t 



dt 2 



Vt, (1) 



where r\ t is Gaussian white noise with zero average, (%) = 0, and correlator {ntVt') = 5(t — t'). 
This is the original Langevin equation with no damping force. The joint random variables 
(xt,vt) evolve in time as 

dv t dx t 

IF = ' 
with initial conditions (xq,vq). Therefore 



dt Vt > -*=*> 



vt = vo + Wt, x t = x + v t + du W u , 

Jo 

where Wt is the integral of the noise, or Brownian motion, hence the process Xt is usually 
referred to as the integral of Brownian motion. 

The statistics of the times of first passage by the origin, and of related quantities, for 
a particle obeying (pi), has been the subject of a long series of works, and is by now well 
understood my- 

More recently, a number of studies have been devoted to survival problems for a ran- 
domly accelerated particle, with particular choices of the boundary conditions at the origin, 
motivated by situations of physical interest. 
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In the partial survival problem, the particle is absorbed with probability 1 — p when crossing 
the origin, and survives with probability p . One is interested in the survival probability, 
that is the probability for a particle obeying (nf) not to be absorbed at the origin up to time t. 
This reads 

oo 

G(i,p) = 5> n p„(*), 

n=0 

where p n (t) is the probability of the occurrence of n zero crossings between and t. At large 
times, the survival probability decays as 

G(t,p) « A(p)t~*<*\ (2) 

t — *oo 

defining the persistence exponent <fi(j>) |6|] and the amplitude A(p). Very recently Burkhardt 
obtained the exact result 



13 p 



4>(p) = - - — arcsin-, i.e., p = 2sin -(1-40) . (3) 



4 2f 2 



n 



6 



He also obtained the persistence exponent O(r) for the inelastic collapse problem, where the 
particle rebounds inelastically at the origin, with coefficient of restitution r < 1, i.e., the 
velocities Vi just before and Vf just after the collision obey Vf = —rvi ||. The result is 



TV 



_(l-40(r))j. (4) 



r 2e M=2sin 
Hence the two exponents are related by 

9(r) = 0(>M), 

as conjectured in ref. [||. 

The purpose of this work is to give new and simple derivations of these results, relying 
on the statistics of the successive times of passage by the origin of a randomly accelerated 
particle. We shall also give an analytical prediction for the amplitude A(p) of the power- 
law (0) [nol. We use the framework of ideas and results contained in the classical work of 
McKean [PJ, complemented by subsequent studies by Lachal M. 

We begin by the partial survival problem. Consider a particle starting from the origin 
with initial speed vq. By definition, 

Pn(t) =V{t n < t<t n+1 ), 

where t n is the time of occurrence of the n-th zero crossing, with probability density f tn (t) . 
Hence, in Laplace space, 



Pn(s) 

where || 



ft As) - f tn+1 (s) 



« ,. 2 f°° 7 sinh(7r 7 /2) K tl (v V8a) 

Jtn i S ) = ~ / "7 ■ 



7T J (2 cosh(7T7/3)) v V8s 

in which K^ is the modified Bessel function. After summation over n > 0, we find 

A, , ^(, n ^2 [°°i 7 sinh(7r 7 /2) K t7 (v Q V8~ S ) \ 

G(s,p) = - 1- (1-J?)- / dj- — — ' \. 5) 

s\ -k J 2 cosh(7r 7 /3) - p vq\/8s I 
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We now have to extract the singular behavior of this expression as s — > 0, corresponding 
to the decay of G(t,p) at large t. To this end we use the integral representation 

THY / 

Ki 1 (x) sinh — = / du sm(-fu) sin(x sinh u) (x > 0), 
2 Jo 

and the integral 

7sin(7u) 3 d sinh (2(1 — <f>)u) 



d/y ■ 

n 2 cosh(7T7/3) — p -v/4 — p 2 du sinh3w 

where <j) and p are related through (y). After an integration by parts, (|[) becomes 

G(s v) MzM 

where 



Jo 



sinh (2(1-^) / n vryl^^ 

-^- cos (t^sinhu) , 1(0) = -g^T^- 



I(s) = / dw coshu cos (i>ov8ssinhu) , 7(0) 



The behavior of the difference 1(0) — I(s) for s — ► is dominated by large values of u. At 
exponential order, setting y = exp(u), we obtain 

1(0) ~m ^Ifdyy- 2 *- 1 [l-cos(tW2^)] = ^ ^ C0S7 ^ ( 2 ^ , 
leading to the result flTj] 

G(i ' p) t^oo 2 ^r (0 + 1) -JT=? [ft) ■ (6) 

For p = 0, this expression reads 



3 ^* 



o(*,o) -«,(*).«_ ^^ )V2t 



£ — >oo 



The persistence exponent of a free randomly accelerated particle, cf> = | P,0l ( see a ls° @)i 
is thus recovered. For p = 1, the exponent = 0, and the right side of (|6|) equals unity, in 
agreement with the identity G(t, 1) = 1. 

The result obtained, equation (|(j), allows the determination of the successive cumulants 
of the number N t of zero crossings between times and t. Indeed, setting p — exp(z), the 
generating function of these cumulants is given by \nG(t,p) = In (cxp(z N t )) . Expanding (||) 
as a power series in z, we thus obtain, 7 E being Euler's constant, 

W) « #(^A-7,)+i (N?)-(N t f « ^_(ln A -7.)+^. ( 7 ) 

t-»oo 2-7T \ 2«q / 6 t^oo 37r y 2«q / 72 

and so on. The mean density of zero crossings in logarithmic time is thus equal to \/3/(27r) H. 
We now turn to the inelastic collapse problem. We consider again a particle starting from 
the origin with initial speed Vq. 
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Following ref. S , we note that equation (pi) is invariant under the change of scale: x — > 
r~ 3 x, t — > r~ 2 t, which implies the rescaling v — ^ r~ x v of the velocity. Hence, performing 
this change of scale immediately after each collision compensates the reduction of speed due 
to inelastic reflection. One can therefore express the elapsed time t\ in the inelastic problem 
in terms of the elapsed time t of a randomly accelerated particle, elastically reflected when 
reaching the origin, as || 



ti = ti+r 2 (t 2 -t 1 )+r 4 {t 3 -t 2 ) + ---+r 2Nt {t-t Nt )= / dur 



2N„ 



where Nt, the number of crossings between and t, and t n , the time of occurrence of the 
n-th crossing, pertain to the elastic problem. According to ref. H, the sum in the right side 



converges if r < 
the finite time 



= n -n-/V3 . 



0.163, demonstrating the inelastic collapse of the particle after 



c = lim ti = h + r 2 (t 2 - h) + r 4 (t 3 - t 2 ) + • • • = / 

t^oo J Q 



dur 2N - 



(8) 



As first observed in numerical simulations || , the tail probability of the collapse time c decays 
asymptotically as 

V(c>t) ~ r~ e W, (9) 

T — >00 

defining the persistence exponent 0(r), the analytical expression of which is given by (0) H. 

We now present new derivations of these results. 

As shown by McKean [nj, a simple consequence of the scaling properties of Brownian 
motion, and of the Markovian character of the couple (x t , v t ), is the existence of the identities 
(in distribution) 

2 

tn = t n — i + u n _ 1 a, 
U n = U n -1 b, 

where u n — \vt n \, a = ti/v 2 ,, b = Ux/\vo\- Hence 

u n = |«o| bib 2 ...b n , 

tn - t n -l = V 2 (6i6 2 . . . b n -l) «n, 

in which the couples (a\,bi), (a 2 ,b 2 ), etc. are independent, with common density Q| 



•• ' V a 



7rv3a 2 



Integrating upon a yields the marginal distribution |I],||] 

The expression of c given by flq) can now be rewritten as 

c = v 2 (01 + r 2 b\a 2 + r 4 b 2 b\az + ■ ■ ■) 
hence (in distribution) 

2 1 2t2 

c = v a + r c. 



(10) 

(11) 
(12) 



G. De Smedt, C. Godreche and J.M. Luck: Randomly accelerated particle 



This is the fundamental equation for the determination of the distribution of the collapse time 
c. The class of problems represented by eqs. (lit) and (Oh has been investigated extensively 
(for a review see e.g. |1T[). In the present case, simple arguments lead to the following results, 
starting from (12). (For a more rigorous treatment, see J12||.) 

The probability density of the random variable c obeys the integral equation 

/ c (c) = / dadbf a>b (a,b) / duf c {u) 6 (c- v„a- r 2 b 2 u) 



= J dadbf^ia^-^U^y^-j. (13) 

Inserting in the last expression a power law of the form / c (c) ~ c~( 1+e - ) for c — ► oo, the powers 
of c on both sides cancel out, and we are left with the condition 

g(@) = ((r 2 b 2 f) = 1, (14) 

which determines the exponent 0(r). Equation (HG) yields 

' ' = 2cos(7rs/3)' 
hence 

■ 9(9) = 2Bin [|(l-4e)]' (15) 

so that (fUJ) is equivalent to (Q). 

The function <?(0) is convex, i.e., g"(&) > 0, and 5(0) = 1. Therefore ( [14|) leads to an 
acceptable solution for the exponent (i.e., a positive one), as long as 

g'(0) = <ln(r 2 6 2 )) < 0. (16) 

If this condition is fulfilled, then the infinite sum ( Jll| ) converges with probability 1, i.e., the 
collapse time c is finite, which is the condition for inelastic collapse. Equation (n5|) yields 
g'(0) = 2(lnr + w/y/3), so that @ is equivalent to r < r c = e _7r/N/5 . 

Let us remark that, in the derivation of (|14|), we implicitly assumed that the argument of 



f c in (13) was dominated by c. This is justified provided that ((v 2 a) 1 ) = ((ii) e ) be finite, 
which actually holds since 0(r) < i for r > 0, as can be seen from (fL5[). 

Contrary to the partial survival problem (see (|6j)), determining the amplitude of the power- 
law (|J) for the distribution of c is, in general, a very difficult task fll],[l^,[l4|. 

Besides persistence exponents, the limiting distributions of observables such as the occu- 
pation time L du (l + signa; u )/2 of the process (|l|), or the instant t^ of last visit of the particle 
to the origin, before time t, arc also of interest. For the case of a particle leaving the origin 
with zero initial speed, Lachal J15J found 

// \ Q f°° p 2z 

F(x) = hm V ( — < x ) = — -r / dz — K Uzx~i)eTi Vfe. 
t^oo y t J 2tt z j z 

Using the integral representation 

/>oo 

K (x)= / due- xcoshu , 
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a simple computation leads to the density 



/w-^w-Gy^r*- "" 







(2 cosh ii-iijvi 2 + cosh u 



i2 



which diverges as f(x) fa Ax s, with A = 327T a (r(^)) /16 ~ 0.224, as x — > 0, and takes 
the finite value /(l) = v3/(27r), at x = 1. This value coincides with the mean density of zero 
crossings in logarithmic time, as expected. 

Determining the limiting distribution of the occupation time of the process (QJ) remains an 
open problem. 

In conclusion, the approach presented here brings, on the problems of partial survival and 
inelastic collapse of a randomly accelerated particle, a viewpoint complementary to that of 
ref. pi], based on Fokker-Planck equations. 



We are indebted to A.J. Bray for having aroused our interest in this problem, and to T.W. 
Burkhardt for stimulating discussions. 
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